Information theoretic aspects in ponderomotive systems 
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We show the possibility to entangle radiation modes through a simple reflection on a moving 
mirror. The model of an optical cavity having a movable end mirror, and supporting different modes 
is employed. The mechanical motion of the mirror mediates information between the modes leading 
to an effective mode-mode interaction. We characterize the modes' entanglement on the basis of 
recent separability criteria. The effect of the thermal noise associate to the mirror's motion is 
accounted for. Then, we evaluate the performances of such ponderomotive entanglement in possible 
applications like teleportation and telecloning. 
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Q , I. INTRODUCTION 

o . . 

Ponderomotive systems are physical systems where the electromagnetic pressure force gives rise to relevant effects. 
The optomechanical coupling between a movable mirror and a radiation field, is realized in such systems when the 
field is reflected by the moving mirror. This coupling was introduced in the context of quantum limited measurements 
ps) ' [1] and then used in interferometric gravitational- wave detection [2] as well as in atomic force microscope [3]. Since 
^ ' then, a wide literature has been devoted to such a coupling. In particular, it has been shown that it may lead to 
nonclassical states of both the radiation field [4,5], and the motion of the mirror [6]. The interest about ponderomotive 
' systems also relies on the possibility to investigate, with them, the tricky borderline between the quantum and the 
LiJ classical world [7,8]. Moreover, recent technical progresses have made this area experimentally accessible [9,10]. 
I The appearance of quantum effects in ponderomotive systems, paves the way to use them also for quantum informa- 

tion purposes [11]. These require as main ingredient the entanglement [12,13]. Furthermore, information processing, 
in the quantum optical framework, can be implemented when applied to continuous quadratures of electromagnetic 
\ modes [14]. Then, the use of a ponderomotive meter for continuous variable entanglement purification has been 
^ ^ investigated in Ref. [15]. Furthermore, the possibility to obtain quantum correlated quadratures of the field refiected 
\ \ by a movable mirror has been proposed in Refs. [16,17]. Here, following the line sketched in Ref. [16], we study 
^ i a ponderomotive system, namely a radiation field refiected by an oscillating mirror, from the quantum information 
d ' perspective. In particular, in Section II, we shall show that the mechanical motion of the mirror mediates information 
^ 1 between the field modes leading to an effective mode-mode interaction. In Section III, we shall characterize the modes' 
C7^' entanglement on the basis of recent criteria [8,18,19]. We shall also account for the effect of the thermal noise associate 
^ , to the mirror's motion. Then, in Section IV, we shall evaluate the performances of such ponderomotive entanglement 
in possible applications like teleportation [20,21] and telecloning [22]. Finally, Section V is for concluding remarks. 



II. A PONDEROMOTIVE SYSTEM 



The model we are going to consider is schematically depicted in Fig.l. It consists of a linear cavity, with an 
oscillating end mirror, plunged in a thermal reservoir at the equilibrium temperature T. This completely refiecting 
mirror, with mass m, can move back and forth along the cavity axes. When the cavity is empty the moving mirror 
undergoes harmonic oscillations at frequency LO„n damped at rate 7„i by the coupling to the external bath. In 
presence of a radiation field, the cavity length varies under the action of the radiation pressure force, which causes 
the instantaneous displacement of the mirror. 
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FIG. 1. A possible scheme implementing the studied ponderomotive system. 
The resonant frequencies of the cavity are calculated at the equilibrium position of the oscillating mirror, resulting 
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where n is an arbitrary integer number corresponding to the index n, c is the speed of light, and L is the equilibrium 
cavity length. We consider the possibility to have several input fields at frequencies won ~ driving the corre- 
sponding cavity modes. In the adiabatic limit in which the mirror frequency is much smaller than the cavity free 
spectral range c/(2L) we can focus only on the driven cavity modes, obtaining the following Hamiltonian 
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where the sum must be intended over the driven modes. Hfree is the Hamiltonian for the free motion of the mechanical 
oscillator (moving mirror) having position Q and momentum P, and of the cavity modes characterized by the ladder 
operators , ^4^ . Instead, Hj^^ive describes the input fields, with amplitudes a™, entering the cavity through the fixed 
mirror whose partial transmission determines the input-output rate 7c. Finally, Hint represents the ponderomotive 
interaction between the mirror and the radiation fields [23]. Such interaction is generated by the radiation pressure 
induced variation of the cavity length, which corresponds to a variation of the frequencies (energy levels) through 
Eq.(l), that is 

ScJcn = gj^ dL = J^Q^ (3) 

with Q = SL <^ L. Since we shall consider few modes whose h ^ c/{2L) differ not too much each other, we can set 
{ucn/L) ~ G, Vn, as the optomechanical coupling constant. 

By using Eqs.(2), and accounting for the losses and the noises, we can describe the complete dynamics of the system 
through the following quantum Langevin equations 

P{t) 



Pit) 



mLolQit) + nGyAi(i)Ar,{i) - 27„P(t) - ^(i) , 



An{t) = -i [lOc 



won) A„(i) + iGAn{t)Q{t) + VT^a™ - ^ A^it) + V7^a™(t) 
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where we have used the replacements An{t) A„e~*'^°"*. Furthermore, a™ are the vacuum noise operators associated 
to the input radiation fields, while ^(t) is the noise operator for the quantum Brownian motion of the mirror. The 
noise correlations are [24,25] 
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with 



where fc^ is the Boltzmann constant. It is worth noting that Eq.(5b) gives the exact thermal noise correlations at 
any temperature T [25]. 

We are now going to study the dynamics of the small fluctuations around the steady state, i.e. the dynamics of the 
operators 

q{t) = Q{t) - X, (6a) 

p{t)=P{t)-y, (6b) 

a„(i) = An - an, (6c) 

where the (classical) stationary values are given by 

X = {Q)ss = I ttn 1^ , (7a) 

y = {P)ss = 0, (7b) 

O^n = (^n)ss = , — / 1 r~r T) (^c) 
UJOn - Wen + Gx 

7c 

the (dimcnsionlcss) overall detuning due to the frequency mismatch and to the radiation phase shift caused by the 
stationary displacement x of the mirror. 

For the sake of simplicity we assume, from now on, symmetric conditions for the various radiation modes, that is, 
A„ = A and a„ = a G 3?, Vn. Then, it is easily recognizable in Eq.(7c) the nonlinear relation between input and 
intracavity intensity field which give rise to the bistable behavior of the system [26]. 

Linearizing Eqs.(4) around the steady state (7) we obtain 

m^^, (9a) 
m 

p{t)= -mul^q{t) + hGY, [a*a„(t) + aai(t)] - 27„p(t) - ^{t) , (9b) 

n 

dn{t) = (^zA - lcan{t) + iGaq{t) + V7^ajf(i). (9c) 

Going into the frequency domain, and eliminating the mirror's variables we are left with a set of 2A'' linear equations 
(A'^ being the number of driven modes, so that n = 1, • • • , N) for the modes quadratures 

Xniu.) = ""^"^ V"^""^ - (IN 

v2 

. (10b) 

v2 

Such equations can be written in compact form as 

iwv(w) = M{lj)v{lj) + ^v'"(w) + s{uj)^(lj), (11) 

where we have introduced the 2N dimensional vectors 

v(l^) - (Xi(c^), yi(w),... , XNiiu),YN{u;)f , (12a) 
v™(c.) = ( Xriu) , Yriuj) , . . . , X^"(c.) , iTM f , (12b) 
s(w) = V2Gx(w) (0, -a,...0, - af , (12c) 

with 

Xi^) = —r^ ^ ' (13) 

' ' m (w^ - a;2 + 2 1 7„ w) 
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the mirror's mechanical response function. Furthermore, A4{(j) is a 2N x 2^^^ matrix written as 
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where Md and Mo are 2x2 matrices given by 
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The useful noise correlations for Eq.(ll) come from Eqs.(5), (10) and read 
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Thus Eqs.(ll)-(17) completely describe the dynamics of the small fluctuations of radiation modes. Practically, we can 
see from Eqs.(9) that the mirror mediates information between the radiation modes leading to an effective mode-mode 
interaction as results from Eqs.(ll), (14) and (15). 



III. OUTPUT FIELDS ENTANGLEMENT 



The above discussed mode-mode interaction will presumably lead to entanglement between intracavity modes, 
which, in turn, should be reflected on the fields outgoing the cavity. On the other hand, only these latter become 
really useful. Hence, we are going to characterize their correlations. First of all we notice, by the the input-output 
theory [24], that 



^ out 
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then, we introduce the hermitian output quadrature operators 
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Their correlations are described by the 2N x 2N matrix 



g= J(v°"*(a;) [v°"*(-a;)]'^ + v''"*(-a;) [w'"'\uj)Y ) , 
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+ i/C(-c.)(v»(-..) [v™(a;)]^) [K.{u)f 
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where 



A:(w) = 7c£"^(a;) - T, (21a) 
C{uj) = iuji - M{uj), (21b) 

with I the identity 2N x 2N matrix. However, the matrix Q is not symmetric and, moreover, it concerns quadratures 
with frequency dependent commutator, i.e., 

( [ ('^) > Ry^' (-o;) ]) =i c{uj) , Vn , (22) 

with c a real positive definite function of frequency cj. Therefore, to easily apply quantum information arguments to 
our system, we construct from Eqs.(20), (22) a symmetric correlation matrix concerning quadratures with canonical 
commutation relations, that is 

^M = ''''^%f'^''^ - (23) 

By virtue of the linearization procedure adopted in Sec. II we have, for each frequency, a multivariate Gaussian state 
completely characterized by Eq.(23). 



A. Bipartite entanglement 

We now restrict the attention to only two mode (A'' = 2) in order to study the entanglement of a bipartite system. 
Practically we consider 



and we introduce the matrices 



Then, Eq.(24) leads to 



V=(c^T^), (26) 



where A and C arc 2x2 matrices. In this case the Simon's criterion [19] is necessary and sufficient for entanglement, 
and, according to Eq.(19) of Ref. [19] we can define a marker of entanglement as 

E = l + {dotAf + Q - |detC|^ - tr {AJCJAJC^j} - ^detA , (27) 

so that, if it goes below 1, the state is entangled. Instead, the product criterion introduced in Ref. [8], and reminiscent 
of nonlocality criterion [27], gives 

E = 4{An+Cn){A22-C22). (28) 

Finally, the sum criterion, expressed by Eq.(3) of Ref. [18], can be rewritten as 

E = ti{A} + tr {Cn} . (29) 

Then, in Fig.2 we report the marker of entanglement E versus u), for the three criteria a) Simon, b) product, c) sum. 
The parameters values arc taken similar to those of the experimental set up of Ref. [10]. They are written in Tab. I. 
Fig.2 shows that in case of no detuning the product and the sum criteria do not reveal any entanglement, while the 
Simon criterion do. This proves the weakness of the the entanglement coming in this case from the interaction of only 
amplitude quadratures as can be evicted from Eqs.(24) and (15). Such type of entanglement, although resistant to 
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thermal effects, is practically useless in information processing like teleportation whose performance are much more 
related to the product and sum criteria [28]. 
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FIG. 2. The marker of entanglement E is plotted versus uj with respect to the three criteria a) Simon, b) product, c) sum. 
Here A = and the other parameter values are given in Tab. I. The dashed lines indicate the limiting value below which 
entanglement is recognized. The three plots remain unalterated in the temperature range T = -^ 300 K. 



TABLE I. Parameter values 
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Then in Figs. 3 and 4 we have only considered the sum criterion and we have shown the beneficial effect of the 
detuning on the entanglement. As matter of fact it allows interaction of both amplitude and phase quadratures as 
can be seen in Eq.(24). In particular Fig. 4 exhibits the presence of entanglement at low frequencies (according to Ref. 
[16]) as well as near the mechanical resonance. Neverthless, the latter turns out to be more sensible to the thermal 
noise. This behavior resembles that of the hght squeezing studied in Refs. [4,5]. 
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FIG. 3. The entanglement indicator E of Eq.(29) is plotted versus u). Here A = —0.1, and solit lines are for T = 0, T = 10, 
T = 50 K from botton to top. The other parameter values are given in Tab. I. The dashed line indicates the limiting value 
below which entanglement is recognized. 
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FIG. 4. The entanglement indicator E of Eq.(29) is plotted versus lj. Here A = 0.1, and solit lines are for T = 0, T = 10, 
T = 50, T — 100 K from botton to top. The other parameter values are given in Tab. I. The dashed line indicates the limiting 
value below which entanglement is recognized. 



B. Tripartite entanglement 

Characterization of multipartite entanglement {N > 2) is a more complex issue [14]. In general, multi-party 
inseparability criteria cannot be formulated in compact form as for the two-party. Here, we consider N = 3, thus the 
matrix 
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Md Mo Mo 

M = \ Mo Md Mo \ , (30) 
Mo Mo Md 

which leads to 

f A C C \ 

V = \ A C . (31) 
A J 

Although for three- mode Gaussian states there exist a necessary and sufficient separability criterion [29] , its violation 
does not necessarily witness genuine tripartite entanglement. However, from the symmetry of matrix (31) we easily 
deduce that the conditions for bipartite entanglement also give tripartite entanglement. As matter of fact a tripartite 
fully inseparable state is that which cannot be separate for any grouping of the parties [29]. But, due to the symmetry 
among the parties, if two of them show entanglement according to Sec. Ill A, then any two of them show entanglement, 
thus revealing a fully inseparable state. 

IV. APPLICATIONS TO REMOTE STATE TRANSFER 

This ponderomotive entanglement find possible applications in quantum information processings with continuous 
variables [14]. Here, we deal with the possibility of using it for remote state transfer. By referring to Fig.l, the modes 
outgoing the cavity can be separated and one of them can reach a sending station while the others reach receiving 
stations. Then, all these modes constitute the quantum channel to exploit for transferring a quantum state from the 
sending station to the receiving ones. We will analyze in detail the case for N = 2, i.e., teleportation [20,21], and 

= 3, i.e., telecloning [22]. 

A. Teleportation 

The standard teleportation protocol for continuous variable [20,21] can be described by a convolution of the Wigner 
functions [30] 

WriP) = J (f^ WsiO K{(3 - , (32a) 
K{P -0=1 Wit -CP- a, (32b) 

where Wr is the Wigner function of the received state, Wg that of the imknown state to be transferred (sent), and 
W that describing the quantum channel between the two stations, i.e., the two entangled modes characterized by the 
correlation matrix (26). Here, small greek letters are for complex variables. 

By Fourier transforming Eqs.(32), we get a simple relation for the characteristic functions namely 



where 



$,(A) = $,(A)^(A), (33) 

K(X) = J d^KK(K) exp(—iKiXi — iK2X2) 

= J d^Kd^/irf^z $(z) exp(— iKiAi — iK2X2) 

X exp{iz • (/ii,-/i2,Ki - A*i,K2 - At2)} , (34) 



where the variables with the subscript 1 (2) represent the real (imaginary) part of the corresponding complex variables, 
and z is a four dimensional real variables vector. Moreover, 



$(z) = exp<; -iz Vz^ } , (35) 
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is the characteristic function describing the two-mode channel, thus characterized by the matrix V given in Eq.(26). 
We also consider a Gaussian state to be transferred, so that 



exp 



-i (Ai,A2) V {XuMf 



(36) 



with V the 2x2 correlation matrix. 

Finally, the fidelity of the protocol, resulting from the overlap between the "r" and the "s" Wigner functions, can 
be written in terms of characteristic functions as 



(37) 



Then, by using Eqs.(33)-(36), we arrive at (see also [31]) 



F = 



1 



4tt 
X exp 



d A exp 



1 



--(Ai,A2)I5(Ai,A2)^ 



— ^ (Ai, — A2, Ai, A2) V (Ai, — A2, Ai, A2) 



^det (21? + W^ATZ + 7^^C + C^TZ + A) 



(38) 



In Fig. 5 we show the teleportation fidelity as function of w. As a state to be teleported we have chosen the coherent 
state for which V ~ diag(l/2, 1/2). In such a case the upper bound for the fidelity achievable with only classical 
means and no quantum resources is 1/2 [28]. Then, we see that this bound is overcame just in correspondence of the 
minima of Fig. 3. Also the behavior of the fidelity in terms of thermal noise reflects that of the entanglement in Fig. 4. 
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FIG. 5. Teleportation fidelity versus lu. Curves from top to bottom are for T = 0, T = 10, T = 50, T = 100 K. Here A = 0.1 
and the values of other parameters listed in Tab. I. The dashed line indicates the classical upper bound. 

B. Telecloning 

As a simple extension of the arguments used for teleportation, we can write the Wigner function of the received 
state (at the two stations) by the convolution [22] 

WriP,v) = J d^^WsiOKiP-C^T^-O, (39a) 

K{f3-^,f^-0= j d^i'WiC -C , ,ri~a, (39b) 

where W is the Wigner function describing the quantum channel between sending and receiving stations, i.e. the 
three entangled modes characterized by the correlation matrix (31). 
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The state at one receiving station can be obtained by tracing the received state (39a) over the other receiving 
station. Due to the symmetry, the two possible states coming out coincide. Thus, we can assume 



Wr{[3) = J d^vWrif3,v) = J <f^Ws{OK{P-0, (40) 

where now 

KiP '0^1 d'v K{P - ^, 77 - • (41) 
By again Fourier transforming Eq.(40), we end up with a relation for the characteristic functions identical to Eq.(33), 

$,(A) = $,(A)i^(A), (42) 

where now 



K{X) = J cPk exp (— zKiAi — zk2A2) 

X j (fiJ.(fCd^z^{z) exp {iz ■ (Aii,-/Z2,Ki - Mi, K2 - /X2, Ci - Mi, C2 - ^2)^} , 



and 



<J>(z) = exp 



4z Vz^ 

4 



(43) 



(44) 



is the characteristic function describing the thrcc-modc channel, thus characterized by the matrix V given in Eq.(31), 
with z a 6 dimensional real variables vector. 

We again consider a Gaussian state to be transferred, as in Eq.(36). Then, the fidelity, being expressed by Eq.(37), 
results, by means of Eqs. (36), (42), (43), (44) as 



F = — I (fx exp 

4 TT 



1 



(Ai,A2)I?(Ai,A2)^ 



X exp 



1 



(Ai, -A2, Ai, A2, 0, 0) V (Ai, -A2, Ai, A2, 0, 0)^ 



^det {2V + -RTAn + TZ'^C + C^U + A) 



(45) 



It practically coincides with Eq.(38). However, in this case, F is limited above by 2/3 [32], due to the no-cloning 
theorem [33]. 

In Fig. 6 we show the telecloning fidelity as function of oj. As state to be telecloned we have choosen again a coherent 
state for which T) = diag(l/2, 1/2). Also in this case the upper classical bound for the fidelity is 1/2 [32]. Then, we 
see that this bound is overcome again in correspondence of the minima of Fig. 4. 




(0/(0m 

FIG. 6. Telecloning fidelity versus uj. Curves from top to bottom are for T = 0, T = 10, T = 50, T = 100 K. Here A = 0.1 
and the values of other parameters listed in Tab. I. The dashed line indicates the classical upper bound. 
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V. CONCLUSIONS 



In conclusion, we have studied ponderomotive entanglement, that is entanglement between radiation modes gener- 
ated by the radiation pressure effects. In doing so wc have also provided a comparison between different entanglement 
criteria. Practically, we have shown that even a classical force, like radiation pressure force, together with macro- 
scopic objects can be used for quantum information purposes. The fidelity achievable in remote state transfer widely 
overcome the classical bounds, thought not reaching the optimal values. Theoretically, this could be obtained by 
an optimization of all involved parameters. However, that would require large numerical resources without adding 
new physics to the problem. So it has been left apart. Instead, we have investigated the role played by the thermal 
noise related to the mechanical motion of the mirror. We have seen that purely quantum effects can survive up to a 
temperature « 10 K. This is within reach in experiments with really macroscopic mirrors [10]. On the other hand, 
the use of micro-opto-mcchanical-systems (MOMS) [34] surely guarantees better performances. 

Beside all that we want to remark that ponderomotive systems have a foundational interest [7] , and, an information 
theoretic approach could help us in understanding the tricky borderline between classical and quantum worlds. 
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